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ABSTRACT 

This thesis presents an analysis of a direct current 
conduction magnetohydrodynamic generator. Laminar, in- 
compressible viscous flow is considered. Experimental 
evidence is offered to show that the laminar assumption 
is valid in the presence of strong magnetic fields. 



Efficiencies are calculated for generators of various 
length and aspect ratio, efficiency curves are presented 
and a method of driving the flow is considered. 
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Chapter I 



INTRODUCTION 

1 . 1 Background 

The minimization of sound transmitted from a submarine 
through the sea to the listening devices of a potential 
enemy is a major problem in submarine design. A large 
fraction of the radiated noise is generated by the submar- 
ine's propulsion plant: the main engines and the screws 

and shafting 

Consequently it would by desirable to devise a propul- 
sion system having neither shaft nor screw and a minimum 
quantity of rotating machinery. Doragh* (l) has shown the 
feasibility of propelling a submarine by means of a sea 
water jet accelerated in a magnetohydrodynamic channel. 
This thesis will be a study of the closely related problem 
of providing the electrical energy necessary to drive the 
magnetohydrodynamic flow. 

The pioneer figure in the study of magnetohydrodynamic 
channel flow was Hartmann (2). In 1937 he made the first 
detailed theoretical analysis of the flow of a conducting 
fluid in a high aspect ratio channel under the influence 
of a uniform transverse magnetic field. Hartmann and 
Lazarus (3), in that same year, reported on observations 
made on mercury flow in a transverse magnetic field. To 
date the principal use of magnetohydrodynamic channel flow 
has been as electromagnetic pumps to circulate liquid 
metal coolants in atomic power installations. 

At the present time there is considerable interest in 
the development of magnetohydrodynamic power generators. 



* The bracketed numbers refer to the Bibliography which 
begins on Page 28 
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Most of these proposed generators make use of a conducting 
gas or plasma as the working fluid. This type of magneto- 
hydrodynamic generating system shows promise of supplying 
large quantities of electric power at thermal efficiencies 
unattainable with present equipment (4, 5, 6, 7) • 

1 . 2 The Object and Scope of This Study 

The object of this study is the design of a liquid 
metal magnetohydrodynamic direct current conduction gener- 
ator compatible with and utilizing the same magnetic field 
as the magnetohydrodynamic pump designed by Doragh. As 
will be shown, the use of a liquid metal, instead of a con- 
ducting gas or plasma, as the working fluid permits the 
design of a complete energy conversion cycle devoid of mech- 
anical pumps and compressors. Since the basic goal of this 
design is to achieve a system free of these mechanical noise 
sources, gaseous working fluids will not be considered. 

The channel under consideration is rectangular in cross 
section and of high aspect ratio. The applied transverse 
magnetic field is uniform and parallel to the shorter sides 
of the channel cross section. Furthermore, the liquid 
metal is assumed to be isotropic, homogeneous and incom- 
pressible 

The study will be restricted to laminar viscous flow, 
since there is little or no experimental data available 
concerning turbulent flow of an electrically conducting 
fluid in a strong transverse magnetic field. In fact, the 
assumption of laminar flow may be quite valid since there 
is some experimental evidence that indicates a transverse 
magnetic field delays the transition from laminar to turb- 
ulent flow. Murgatroyd (8) in his experiments with mercury 
found that laminar flow was stable when the hydraulic 
Reynolds number was less than 225 times the Hartmann number. 
Lock (9) hypothesized that laminar flow would be stable for 
a Reynolds number less than 42,000 times the Hartmann number 



3 

for larger Hartmann numbers. Since this study will be deal 
ing only with large Hartmann numbers, both theory and exper 
irnent predict laminar flow. 
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Chapter II 
PROCEDURE 
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2.1 Basic Equations 

The fluid under consideration is homogeneous, isotropic, 
non-magnetizable , non-polarizable and incompressible. It is 
characterized by a permeability p 0 , an electrical conductivity 
a mass density p and a dynamic viscosity • The equations 
that govern the motion of such a fluid are the magnetohydro- 
dynamic Navier-Stokes equation. 

= -grad p + v^v + JXB, (2.1.1) 

the incompressibility condition 

div v = o, (2.1.2) 

Maxwell's equations 

curl H = J, 
curi E = , 

and 

div H = o; (2.1.5) 

and the constitutive relation 

J = CT(e + v X pH). (2.1.6) 

Where H is the magnetic field intensity, J is the conduction 
current density, E is the electric field intensity, v is the 
fluid velocity and P is the pressure. Rationalized MKS units 
are used in all equations and all equations are written with 
respect to a fixed coordinate system 



( 2 . 1 . 3 ) 

(2.14) 



5 



2. 2 Pressure and Velocity Relations 

The geometry of the generator to be considered is il- 
lustrated in Figure I. Fluid flows in the z -direction 
with a velocity that is a function of both x and y. A 
uniform magnetic field H is applied in the x direction. 
The channel walls at y = - w/2 are highly conducting 
buses that are maintained at a potential difference by 
the external constraint imposed by the load resistance. 
The walls at x = - d/2 are non-conducting. The electric 
field E resulting from the transverse potential gradient 
is assumed to be in the negative y direction, and the 
current density J is assumed to be in the y direction. 

The channel cross section is shown in Figure II. 

Load current I is assumed to flow symmetrically in the 
external circuit This insures that the magnetic field 
resulting from the flow of load current will be in the 
z direction and so will not affect the generator oper- 
ation. End effects will not be considered at this time. 

Under these assumptions , the z component of Equation 
(2.1.1) becomes 



0 = - - - J y B x 



bz 



tr 



( 21 . 1 ) 



The elimination of J by use of Equation (2.1.6) yields 
-lii + Vv - f[_B 2 v + £_( EB -1 !_!) = 0 * (2.2.2) 

iy 2 l ^ ^ 

This equation can be written as 



* From this point on indices indicating vector directions will 
be omitted. All directions are as defined in Figure I. 
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Channel Cross Section and External Electric Circuit 
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* 2y , > 2 v 

c>x 2 b y 2 



k 2 v 



+ P = 0 



where 

k 2 = B 2 (d/2) 2 
-(d/2) 2 




( 2 . 2 . 3 ) 



(2. 2.14-) 



M = B(d/2)^JiL 



(Hartmann number) 



and 



? = (EB - — p ) 



(2.2.5) 



( 2 . 2 . 6 ) 



A solution to Equation (2.2.3) of the form 

V =21 v ?n - 1<X) ° 0S ( 2n ~ l)lty dy (2.2 7) 
<v*- r c w 



is desired, where 



2n 



_! ( x ) = £ 
w 



^L 

V V cos 

-w„ 



t 2 "- 1 ^ dy 



( 2 . 2 . 8 ) 



w 



If Equation (2.2.3) is multiplied by the function 
cos (2n-l )ny an( j resulting equation integrated 



w 



over the Interval y = -w/2 to y = w/2, and the boundary 
condition, v(x, -w/2) = 0 is applied, the result is 
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b V Pn- 



2n-l - 



M‘ 



(<a/2)‘ 



t~ (2n-l )tc 1 2 v 2n _ + 

+ L w J (2n-l)rc 



(2n-l )i 

(2.2.9) 



or 



'b v 



2n-l 






n + 1 



= 0 



^x 2 



(2n-l)n 



( 2 . 2 . 10 ) 



where 



M 2 = M 2 0 + ( 2n " 1 ) ,lc 

n (d/2) 2 w 



( 2 . 2 . 11 ) 



Define the function 

6 2n-l " "nSn-l + 



n + 1 



(2n-l )it 



( 2 . 2 . 22 ) 



then Equation (2.2.10) can be expressed as 



2 

^ ^ 2n-l - Mn 2 £ 2n _ ± = 



0 



i - 2 



X 



(2.2.13) 



This has the solution 

6 2n _ 1 = Acosh M n x + B sinh M n x (2.2.1k) 

where A and B are arbitrary constants to be evaluated from 
the boundary conditions. 

Because of the symmetry of the problem, v, and hence 
6 2n-l must be an even function of x. Therefore B = 0. 

A may be evaluated from the requirement that v(-d/2,y) = 0. 
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The solution for 6„ . is then 

2n-l 

_ -L|.P (-1 ) n + 1 
°2n-l 

(2n-l)it cosh M n^ 

2 



(2.2.15) 



and 



4P(-l) n + 1 n cosh Mx 
2n-l 1 ~ S— 



(2n-l)itM 



n 



cosh M d 
n 



(2 2.16) 



whence 






V = 



11 




p '=» 



(- 1 ) 



n + 1 



(2n-l)M n ‘ 



, cosh M x 
1- n 



cosh M d 
n 



The expression for the mean velocity _v can be found by 
integrating Equation (2.2.17) across the channel with the 
result 



v 




(2n-l )^M 



2 W 2 



n 



tanh 

M d 
n 




(2. 2.18) 



Equations (2.2.6) and (2.2.18) can now be combined to 
give the pressure gradient down the channel. 



